CO 


NASA TECHNICAL MEMORANDUM NASA-TM-76658 19820013356 


ifySPi TlTi - 'y la 

NASA TM-76558 


76 00168 9026 


THE INFLUENCE OF A LOCAL WALL DEFORMATION ON THE 
DEVELOPMENT OF NATURAL INSTABILITIES IN A LAMINAR 

BOUNDARY LAYER 


S. Burnel, P. Gougat, and F. Martin 


Translation of "Influence d'une deformation locale de paroi sur le 
developpement des instabilites naturelles d'une couche limite 1am- 
inaire", Journal de Mecanique, Vol . 16, No. 3, 1977, pp. 417-439 


> 

LIBRARY COPY 

JUN O L 1981 


LANGLEY RESEARCH CENTER 
LIBRARY, NASA 
HAMPTON, VIRGINIA 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 JUNE 1981 


stampard title pace 


)• Report No. 2. Oovornmertt Acc.taion No. 

NASA TM-76558 

3’, PscipientsCo*oiogNo. * 

. ^ 

4. Title ond Svbfifle 

THE INFLUENCE OF A LOCAL WALL DeFOR!4ATIO^ 
ON THE DEVELOPMENT OF NATURAL INSTABILI- 
TIES IN A LAMINAR BOUNDARY LAYER 

5. Report Dote 

June 1981 

6. Performing Orgonizotion Code 

7. AutKor(t) 

S. Burnel (University of Orleans-La 
Source) and P. Gougat and F. Martin 
(CNRS Aerothermics Laboratory) 

6. Performing Orgonizotion Report No. 

10. Work Unit No. 

\ 

9. Performing Orgonizotion Nome ond Address; 

Leo- Kanner Associates 
Redwood City, California 94063 

11. Controct or Gront No. 

NASw-3199 

13. Type of Report ono Period Covered 

Translation 

13. Sponsoring Agency Nome ond Address 

National Aeronautics and Space Adminis- 
tration, Washington, D.C. 20546 

14. Sponsoring Agency Code 

15. Supplementory Notes 


Translation of "Influence d.'une deformation locale de pafoi 
sur le developpement- des instabilites naturelles d'une couche 
limite laminaire"’, Journal de Mecanique, Vol. 16, No. 3, 1977, 
pp. 417-'439. 


16. Abstroct 


Abstract.— T he natural instabilities which are propagating in the laminar boundary 
layer of a flat plate are composed of intermittent wave trains. 

A spectral analysis of these fluctuations permits to determine their frequency range 
where always arises a frequency / and the harmonic 2 / only if there is a wall deformation. 
This analysis gives also the amplitude modulation spectrum of the instabilities. Plots 
of the evolution of their power spectral density are compared with the numerical results 
got from the resolve of the Orr-Sommerfeld equation, while the harmonic is related 
to a micro-recirculaling flow’ near the w'all deformation. 


17, K«y Wordt (Selected by Aulhor(s)) 


\ 


18. Oisfribution Srot«m«nt 

Unclassified-Unlimited 


19. Security Clossif. (of this report) 

20. Security Cloieif. (of thi i poge) 

21. No. of Poges 

22. 

Unclassified 

Unclassified 




nasa-hq 


y-fY) - 76? 3 f> J 


THE INFLUENCE OF A LOCAL WALL DEFORMATION ON THE 
DEVELOPMENT OF NATURAL INSTABILITIES IN A LAMINAR 

BOUNDARY LAYER 

S. Burnel*, P. Gougat**, and F. Martin 

Abstract /417 

The natural instabilities propagated in the laminar boundary layer 
of a flat plate are composed of intermittent wave trains. 

Spectral analysis of these fluctuations can determine their 
frequency. There always arises a frequency f, but the harmonic 2f 
only appears if there is a wall deformation. Spectral analysis also 
gives the amplitude modulation spectrim of the instabilities. Plots 
of the evolution of their power-density spectrum are compared with the 
numerical results obtained from solution of the Orr-Sommerfeld 
equation, while the harmonic is related to a microflow along the 
wall deformation. 

Introduction 

In the boundary layer of a flat plate, the transition from a 
laminar state to a turbulent one starts with the development of 
unstable frequencies. This harbinger of transition has been studied 
by Tollmien [1] and Schlichting [2] : By introducing two-dimensional /418 
velocity perturbations, defined by a current function, into the 
Navier-Stokes equations, the Orr-Sommerfeld equation is obtained. 

Its solution supplies the elements for the theory of stability: 

Depending on the magnitude of its reduced frequency w* and the 
Reynold's munber of the flow the perturbation introduced into 

the boundary layer is either amplified or attenuated by a coef- 
ficient a^. The linear amplification in the first stages of the 

*University of Orleans-La Source, Applied Mechanics and Energetics 
Laboratory, 45100 Orleans, France 

**CNRS Aerothermics Laboratory, 4 ter, route des Gardes, 92190 
Moudon, France 

^Numbers in the margin indicate pagination in the foreign text. 
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process then gives rise to spots of turbulence whose coefficient of 
intermit fence grows until an actual state of turbulence is reached. 

It seems that the entire transition phenomenon is ruled by 
the development of such instabilities. The first experimental 
verification of this was made by Schubaure and Skramstad [3] , who 
used a vibrating band to introduce sustained perturbations in the 
boundary layer. More recently, Jordinson [4] took into account 
the thickening of the boundary layer and studied the influence of a 
pressure gradient on the development of instabilities. These 
different comparitive studies dealt with artificial initiations of 
the transition phenomenon. While we have constantly used such 
work as a point of reference, our own study in fact concerns the 
influence of a wall deformation on the position of natural transition. 
The complexity of the phenomenon has led us to look at both the case 
of a wall deformation and that of a flat plate at the same time. 

In the course of our study of natural instabilities, we en- 
countered an unexpected problem — the appearance of a harmonic 
of the instabilities at a right angle to the wall deformation. 

2. Experimental Equipment 

2.1. Wind Tunnel 


Measurements were made in a subsonic wind tunnel of the Eiffel 
type. The test stream had a square cross section (0.5 m on a side) 
and a length of 1.5 m. It was preceded by a convergent nozzle with 
a contraction ratio of 9. The speed of the flow could vary from 
4 to 30 m/sec. For speeds of between 10 and 25 m/sec, the relative 
level of longitudinal speed fluctuations. /u'^/U ^ ,, remained 
basically constant and equal to 0.4%. In addition, spectral 
analysis of the fluctuations in the free flow speed revealed a flat 
spectrum without any singularity that could prematurely initiate 
transition. A weaker preturbulence might push transition slightly 
upstream without changing the structure of the phenomenon. Such 
a shift would not interfere at all with the comparative study of 
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the development of instabilities in the presence of a wall singularity. 


An incremental two-directional displacement mechanism made it 
possible to position hot film sensors with 0.1 mm steps in the 
ordinate and 1 mm ones in the abscissa. 

2.2. Deformable Plate 


The study was conducted in the boundary layer of a deformable 
flat plate. To avoid spurious vibrations, the plate was placed on 
four supports held by a chassis independent of the wind tunnel struc- 
ture. It was preceded by an elliptical leading edge with an aspect 
ratio of 10. Studying the effect of a local wall deformation on the 
boundary layer was greatly facilitated by the machined deformable 
plate. The plate served as a cover for a metal frame in which 
three independent, sealed cavities were bored. In the lower section 
of each cavity were two orifices for increasing or reducing the 
air pressure to produce concave or convex deformations in the plate. 
The amplitude of the deformations was measured by a capacitive 
displacement pickup. The beginning of the first cavity was located 
at an abscissa of 178 mm from the leading edge. The length A of 
each cavity was 100 mm. The abscissa relative. to the first cavity 
will be referred to as x* = x - x^. Only the first cavity was 
used in the present study, and we worked with constant pressures, 
which led to unchanging wall deformations. The amplitude of the 
wall deformations did not deviate by more than 10% from a sinusoidal 
variation of the type y = a [1-cos (2irx*/X ) ] . The maximum amplitude 
2a will be used to characterize the attained deformation at x* = X/2. 
Finally, note that the two-dimensional character of the deformation 
was conserved over 90% of its span. 

2.3. Hot Film Anemometer 


The measurements of average velocity and of velocity fluctuations 
were made with the help of an anemometric circuit. The sensitive 
element was either a wire or a film. The results presented here were 
obtained with a platimum film depositied on a quartz rod 25 pm 
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in diameter. We simultaneously used a split film sensor to find the /420 

direction of the velocity vector in the boundary layer. The sensor 

was made up of a quartz bar 150 ym in diameter and 2 mm long on which 

two platimum films were deposited. The plane dividing the two 

films was parallel to the plane of the plate. The films were controlled 

at constant resistance by two independent aneraometric circuits. In 

two-dimensional flow, the ratio between the heat fluxes and 

released through convection at each film made it possible to measure 

the angle made by the velocity vector with the plane dividing the 

sensor. The sum of the heat fluxes on each film gave the magnitude 

of the velocity vector. 

The sensor was adjusted through the choice of functional 
resistances for each film. These resistances were chosen so as to 
make the ratio ^p/^-Q equal to one for an angle of incidence a of 
zero [5]. Under such conditions, for a = 90°, the ratio 
remained constant for speeds ranging from 5 to 16 m/sec. 

We calibrated the differen- 
tial probe for an angle a of 
10°: It is interesting to note 

that the variation of the fluxes 
and followed a Collis and 
Williams type of power function 
as if only a single cylindrical 
film were involved. In figure 1, 
we have plotted the flux of each 
film as well as the total flux 

Calibration of the /421 

sensor, giving the variation of 
the ratio R between the fluxes 
and il>g as a function of the 
angle of incidence a and the 
Reynolds number, has been accomplished by Olin and Kiland [6] . We 
adopted the following approximate independent quadratic equation 
to analyze the measurements: 



Figure 1 

Exterior Flow Speed Ue vs. 
Total Convection Flux 'J’a+^’b 
and also vs. the Flux of Each 
of the Probes Films ^b) 
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K0-R90 V 90; ’ 

with R„ and representing the values of the ratio R for angles 

u y u 

a of 0 ° and 90°. 

2.4. Frequency Analysis 

The power-density spectrum of velocity fluctuations was 
determined by means of a real time analyzer. 

The apparatus included an intake memory which received the 

signal in the form of a digital sample. The rate the memory was 

reread at was greater than the data entry rate. The signal samples' 

time scale was compressed before the data was fed into a heterodyne 

filter. This allowed the exploration of frequencies to be accelerated 

while at the same time maintaining a high degree of statistical 

precision. As an example, for an analysis ranging from 0 to 

1,000 Hz, the sampling frequency would be 3,000 Hz, and the time 

it would take to fill up the memorv would be 0.2 sec. With a time 

-3 

compression factor of 1,000, in 40 x 10 sec we would have an 
instantaneous analysis of 200 points with an effective resolution of 
9.4 Hz. (This figure takes into account the frequency sweep and 
the weighting function for signal truncation effects.) 

The power-density spectrum was obtained by numerical integration 
of a great number of instantaneous spectra. Our measurements were 
made by accumulating 2,048 spectra, which corresponds to a sampling 
time of 6 min. 

The use of a tape recorder in conjunction with the real time 
analyzer provided us with some interesting possibilities. The 
velocity fluctuation signal was recorded and played back at a slower 
speed before being fed into the anlyzer's memory. We thus had a 
simultaneous oscillocope display of a portion of the instantaneous 
sign.rl and of the instantaneous spectrum of this sample. The 
change in speed brought about by the tape recorder also enabled 
us to watch the signal go by and block in the memory the most 
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significant transient portions. 


3. Influence of the Wall Deformation on the Exterior Velocity Gradient /422 


3.1. Measurement of the Velocity Gradient 

The development of the 



Figure 2 

Polhausen Parameter Curves 
Showing Experimental Values 
Obtained as a Function of 
X and X* for Different Wall 
Deformation Amplitudes 
(Uoo = 16 m/sec) 


is then of the form: 


boundary layer is a function of 
the boundary conditions. It is 
therefore as necessary in flat 
plates as deformed ones to 
measure the longitudinal 
velocity gradient at the edge of 
the free flov; boundary layer. 

The results are presented in 
figure 2, where the gradient is 
characterized by the value of 
the Polhausen parameter A = 

(6^/v) (dU^/dx) . The experimental 
curves intersect at x* = 0 and 
X* = 50 mm. This result is 
explainable when the deformation's 
profile is sinusoidal. The 
potential flow along the profile 


Ue , 2k a _2„v/x 2k X* 
— = 1 + e , 

U, X X 

from v.’hich . .* 

A = Kae~^’"'^*'sin 


which corresponds to a bundle of curves cutting the A = 0 axis at 

the points x* = 0 , x* = A/2, and x* = A . The experimental curve 

bundle differed slightly from this because the deformation of the 

plate had a finite length A. In addition, the elliptical leading /423 

edge induced a velocity gradient which modified the potential flow 

expression . 
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3.2. Transition Position 


Previous experimental work showed that the phenomenon of 
transition on a flat plate is not characterized with enough precision 
by the evolution of the average velocity profile, which gradually 
changes from a Blasius profile to a turbulent one [7]. In contrast, 
the study of how longitudnal fluctuations develop provided a 
"transition criterion" that was a little more precise. 




Figure 3 

Velocity Fluctuation Rate 
vs. Ordinate on a Flat 
Plate for Various Abscissas 
(U^ = 12.5 m/sec) 


Figure 4 

Influence of the Wall Defor- 
mation on Transition Position 

(U = 16 m/sec) 

00 


In figure 3, we show how the /u ' ^/Ug level changes as a function 
of the ordinate for various abscissas. Both the maximum value of 
>^u ' and the ordinate at which it is achieved depend on the 

abscissa. To avoid the difficulties of a three-dimensional 
representation ( / a ' ^ /U^ = f(x,y)), we have characterized each of the 
/u'2/Ug = f(y) curves by the value of its maximum, which we have 
labled 

' e max 


The variation of the quantity (/u'^/U ) as a function of 

0 rri3.x 

the abscissa is presented in figure 4 for both the flat plate and 
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tv70 different wall deformation amplitudes. These curves have in 
turn a maximum whose abscissa could be adopted as the criterion for 
transition . 

Thus, a negative wall deformation with a -1 mm amplitude causes 
transition to shift forward by 37%, Note that there is only a 3% 
difference between the forward shift due to positive deformations 
(+1 mm) and negative ones (-1 mm) . 

In fact, the abscissas at which this criterion is applied 
correspond to states of transition which are already advanced, and 
the overall measurements only illustrate the significant effect a wall 
deformation of small amplitude has on the phenomenon. 

In order to study pretransitional mechanisms, we had to 
investigate the detailed structure of the instantaneous signal in 
the zone where the Reynolds number was close to the critical value. 
This enabled us to discover the frequency range of the natural 
instabilities and their coefficient of amplification. 

4. Observation of the Velocity Fluctuation Signal 

4.1. Elimination of Low Frequencies 

Oscilloscopic observation of the instantaneous signal showed 
the existence of low frequencies (below 200 Hz) . We attempted to 
describe the changes in the effective value of such frequencies. 

They were present from the leading edge onward and existed in the 
laminar zone. Their variation as a function of the abscissa did 
not seem to exhibit any distinguishing feature. It seemed possible 
to explain them by the unstable position of the stop point on the 
leading edge. 

We will limit ourselves to mentioning their existence. In so 
far as the low frequencies did not appear to influence transition, 
the rest of the study used a filtered instantaneous signal. 
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4.2. Integrated Spectra 


Thus, after filtering out the low frequencies, the velocity 
fluctuation signal was fed into the analyzer. The signals' 
power-density spectrum is presented in figure 5 for an abscissa 
X* = 50 mm and an ordinate y = 0.1 mm, which is where the signal to 
noise ratio was the best: 

— Case a refers to a flat plate. The frequency ranges from 
300 to 700 Hz, and, for a Reynolds number of 630, the maximum 
energy level within this range is located at the frequency f^ = 

500 Hz. 


-- Case b refers to a deformed plate. In addition to the 
energy contribution at 500 Hz, a second maximum is located at £2 = 
1,000 Hz for the Reynolds number under consideration. 

4.3. Instantaneous Signals and Their Spectra 


Detailed observation of the instantaneous signals was obtained 
with the aid of the real time analyzer after speed reduction by'means 
of a tape recorder. The instantaneous signal sample was put in the 
memory and then analyzed. Simultaneous displays of the instantaneous 
signals and of their spectra are shown in figure 5: 

-- Case c corresponds to a signal sample obtained on a flat 
plate: The instantaneous signal was a relatively pure sinusoid, 

modulated in amplitude. The instantaneous spectrum, like the 
integrated one, exhibits a maximum energy level at a frequency of 
500 Hz. 


— Cases d, e, and f correspond to three samples involving 
deformed plates. The instantaneous signals were no longer pure 
sinusoids. For cases d and e, the instantaneous spectra simultaneously 
exhibit a frequency f^ of 500 Hz and a frequency £2 of 1,000 Hz. 

In contrast, case f's spectrum only exhibits the 1,000 Hz £2 
frequency . 


9 



Figure 5 

Integrated Spectra of a Portion of the Signal Obtained: 
a) Flat Plate, 2a = 0, at x* = 50 mm and y = 0.1 mm 
b) Deformed Plate, 2a = -1, at x* = 50 mm and y = 0 . 1 mm 
Unstable Bursts Characteristic of the Hot Wire Signals 
and their Instantaneous Spectra (Uoo = 16 m/sec) : 
c) Flat Plate; d,e,f) Deformed Plate 


The structures of instantaneous signals obtained with a 
flat plate and with deformed ones v;ere similar . In both cases, 
the fundamental frequency appeared in the form of intermittent 
wave trains. In deformed plates, the fundamental was altered by 
a harmonic, but its structure remained the same (frequency, inter- 
mittent character) . 

4.4. Identification of the Harmonic 


The signal samples reproduced in figure 5 show that the f 2 


/42" 
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The harmonic of the insta- 


frequencies were twice those of 
bilities was thus involved here. (Note that only the first harmonic 
existed.) This was confirmed for other free flow speeds and other 
abscissas: When the Reynolds number was changed, the f^ frequency 

varied, but the f 2 frequency remained double f^. What is more, 
examination of the instantaneous signals in figure 5 shows that the 
energy present in the spectrum at f 2 is due to a perturbation 
grafted on to the sinusoid of frequency f . The fact that only the 
first harmonic existed eliminates the proposition that the electronic 
signal processing circuits were saturated. For the same amplitude 
wall deformation and for a given sensor position, all the samples 
were not contaminated in the same way by the f 2 frequency. In 
particular, we see in figure 5 that the contribution of f 2 was greater 
in sample e than in sample d. Sometimes the f^^ frequency disappeared 
completely to the benefit of frequency f 2 • This was the case in 
sample f . 

4.5. Identification of the Fundamental Frequency 

Obremski ' s and Morkovin's theoretical results [8] concern solving 
the Orr-Sommerfeld equation so as to trace the equal amplification 
curves in the w*, plane for the case of a spatial variation in 

the instabilities. At a given Reynolds numbei; it is possible to 
calculate the range of unstable frequencies having maximum ampli- 
fication . 

The experimental results concerning the f^ frequency on a flat 
plate for different Reynolds numbers were in good agreement with the 
theoretical ones (figure 6) . However, it is interesting to return 
to the instantaneous signal samples of figure 5 in order to note that 
the period of the phenomenon seems constant although the analysis 
resulted in a wide spectrum around f^. We have tried to interpret 
this phenomenon in terms of amplitude modulation. The details are 
given in figure 7. Let us start with a section of the instantaneous 
signal (figure 7a) and its spectrum (figure 7a') and look for an 
analogy with a simple model of a pure sunsoid of frequency f whose 
amplitude is modulated by a low frequency f (figure 7b). The 


/428 


11 


spectrum shows a line at 
frequency f and two modulation 
lines at f ± f (figure 7b'). 

The conventional approach in signal 
processing for isolating the 
f ' frequency is to carry out a 
full wave rectification of the 
signal (figure 7c) . The spec- 
trum of the rectified signal 
includes a line at a frequency 
of 2f, two modulation lines at 
2f ± 2f, and at all the higher 
harmonics as well, and finally 
a line at f (figure 7c'). The 
model can in fact allow for mod- 
ulation by narrow band noise 
(figures 7d and 7d'), and recti- 
fication would result in iso- 
lating a continuous modulation 
spectrum (figures 7e and e'). 

The same procedure of full wave rectification could be applied to the 
sample of the physical signal (figure If.) . Analysis of the 
rectified signal uncovers a whole range of low frequencies con- 
stituting a continuous modulation spectrum (figure 7f). We can thus 
explain why the instantaneous signal samples v;ere made up of pure 
sinusoids while their spectra only contained a single line. In 
fact, it is difficult at the present state of signal processing to 
decide on the origin of the modulation spectrum: Two phenomena are 

simultaneously at play, the actual modualtion of the unstable bursts /430 
and the intermittence function for these bursts. We have hypothesized 
that these two phenomena are provoked by the low frequencies present 
in the boundary layer, but the correlation between the modulation 
spectrum and the spectrum of low frequencies cannot be documented 
except by modifying the low frequencies' power-density spectrum. 



Figure 6 

Frequency vs. Reynolds Number 
in Relation to Thickness 
of Boundary Layer Displacement 
j — Uoo6 1 /v 

, Theoretical Maximum 

Amplification Curve 
O, Experimental Points 
2a = 0 


Key: a) Theory 

b) Experiment 
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Figure 7 

Amplitude Modulation by a Low Frequency Spectrum 
Time Scale: 8 x 10”^ sec/square for Instantaneous Signals, 

100 Hz/square for Spectra 

a) Unstable Burst on a Flat Plate for a Reynolds Number of 1110 
b) Sinusoid with Amplitude Modulation 
c) Signal b after Full Wave Rectification 
d) Sinusoid Modulated by Narrow Band Noise 
e) Signal d after Full Wave Rectification 
f) Unstable Burst in a after Full Wave Rectification 
a ' , b ' , c ' , d ' , e ' , f': Spectra of Signals a,b,c,d,e,f. Respectively 
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5. Measurement of Energy for the Fundamental Frequency 


5.1. Variation as a Function of the Ordinate 


VJe have now established that for a given Reynolds number, the 
instabilities have a fixed frequency f^ and that this frequency 
corresponds to the power-density spectrum maximum for the velocity 
fluctuations in the range of analysis under consideration. V7e 
will designate the energy contribution made by the fluctuations 
of frequency f^^ as ( /u ' ^/Ug) • 

We will only examine the variation in this quantity at a single 
abscissa (x* = 100 mm) , and only for two deformation amplitudes 
(2a = +1 mm and 2a = -0.5 mm) , but the phenomenon is of the same type 
for a flat plate. 

Case of the Positive Deformation 2a = +1 mm 



Figure 8 

Energy of fi Insta 
bilities vs. y/6 
and Profiles of Cor 
responding Average 
Speeds 

Uco = 16 m/sec, x*=100 
6 2 . 51 , 


mm 




of the first maximum. 


The variation in energy at frequency 
f^ as a function of y/6 is shown in figure 
8. A first maximum exists at y/6 = 

0.18. This maximum can be attributed 
to the presence of a point of inflection 
in the average velocity profile. According 
to Karman-Pohlhausen , the velocity profile 
can be represented by a fourth order poly- 
nomial approximation in y/6 . The poly- 
nomial's coefficients depend on the 
parameter A, characteristic of the exter- 
ior velocity gradient. The point of 
inflection on the velocity profile 
(d^ (U/U^) /d (y/6 ) ^ = 0) is obtained for 
y/6 = -A/(12-2A). For our experimental 
conditions, x* = 100 mm, 2a = +1 mm, one 
finds that A = -3.5 and thus y/6 = 0.184, 
which is little different than the ordinate 
This correlation between the amplification of 
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of the instabilities and the point of inflection on the velocity 
profile was proposed by Tollmien to be a singularity of the Orr- 
Sommerfeld equation. 


The second maximum occurs at y/6 = 0.35, which is in good 
agreement with the experimental results of Ross, Barnes, and 
Burns [9]. Note that not only the value of this maximum, but also 
its ordinate in the boundary layer varies as a function of the 
Reynold's number. However, in all cases the maximum's ordinate 
is such that the average speed at this point is U/Ue = 0.4, which 
corresponds to another singularity in the Orr-Sommerf eld equation, at 
which the velocity is equal to the phase velocity of the instabilities 
Thus, it is the value of the maximum that is useful to us in char- 
acterizing the state of instability amplification at a given abscissa. 

A third maximum exists in the neighborhood of the boundary 
layer's edge. We do not attach any importance to this maximum 
because it concerns an input of energy at the f^^ frequency resulting 
from the general return of the spectrum to the continuous spectrum 
of free flow. 

Case of the Negative Deformation 2a = -0.5 mm 

The results are presented in figure 8. 

There is no point of inflection on the average velocity 
profile and therefore no first maximum at y/6 = 0.2. This is due 
to the fact that for the abscissa and deformation under consideration, 
the A parameter is positive. 

The maximum energy level is located at the point where the 
speed U/U^ = 0.4. 

Thus, the variation of the power-density spectra of the 
fj^ instabilities along the normal to the wall depends on the sign 
of the exterior velocity gradient. However, the variation exhibits 
an extremum characteristic of the state of instability amplification 
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in both cases presented here. We can now justify the choice of 
these two deformation amplitudes for figure 8. We chose a positive 
and a negative amplitude leading to the same state of amplification 
at X* = 100 mm. We thus showed that the existence of a first 
maximum is independant of the state of instability amplification and 
depends only on the presence of a point of inflection on the 
velocity profile. 

5.2. Variation as a Function of the Abscissa 



Figure 9 

Amplification Coefficients 
of f^ Instabilities vs. 
Exterior Velocity Gradient 
(Uoo = 16 m/sec) 

Key: a) Velocity Gradient 


The variation of the maximum 
level of instability as a function of 
the abscissa is traced in figure 9. 
The curves in the figure express the 
spatial amplification of the insta- 
bilities for a flat plate and for two 
different defoinnation amplitudes. 

The amplification curve for the 
flat plate is monotonic: 


Its slope represents the amplifica- 
tion factor. For an abscissa x* = 

50 mm corresponding to a Reynolds 

number N„. = 870, the exoerimental 

Ko 2 

amplification coefficient -a.6,/N„. 

1 J- xxO 2 ^ 

at the reduced frequency oj* = 1.2 x lO"*' 

-5 

equaled 2.1 x 10 and did not differ 
by more than 3% from the the theoretical 
value obtained from the Obremski- 
Morkovin tables . 


The curves for the defomable wall exhibit localized changes 
in curvature at abscissas which depend on the sign of the deformation : 


— For a positive deformation, the favorable gradient (A > 0) 
intervenes at the beginning of the amplification and attenuates the 
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instabilities. Their energy level at x* = 60 mm is less than at 
X* = 30 mm. 

— For a negative deformation, the instabilities are already 
greatly amplified by the time the favorable gradient (A > 0) 
intervenes. The gradient can only reduce the amplif ication very 
slightly and does not succeed in attenuating the instabilities. 

The amplification curves in figure 9 are not traced beyond the 
abscissa x* = 100 because in the presence of a wall deformation of 
either sign, the unstable bursts exhibit a wide spectrum corresponding 
to a complete breakup, and it becomes impossible to characterize 
the evolution of transition by the selective amplification of one 
unstable frequency. 

At a right angle to the deformation in a given section, the 
effect of the exterior velocity gradient depends on the energy 
level attained by the instabilities in the section. There is there- 
fore no direct corrspondence between the state of amplification and 
the local value of the A parameter. 

6 . Energy Measurement for the Harmonic 

6.1. Variation as a Function of the Ordinate 


In figure 10 we have plotted the relationship between y/6 
and the energy level of the fluctuations in relative velocity at a 
frequency f, i.e. the quantity ( /u' 2/U^) ^ for frequencies f^^ (curve 
a) and £2 (curve b) . These values were obtained for a wall de- 
formation amplitude 2a = -1 mm and an abscissa x* = 50 mm. Curve 
b shows that the unstable frequency £2 had maximum energy at a 
reduced ordinate y/6 of 0.1, which corresponds to a small distance 
from the wall (between 0.1 and 0.2 mm) . Curve c represents the 
variation as a function of y/6 of the ratio ( /u ' 2 /U ) /( /Q ' 2 /u 

^ dQ G f ^ 

which was at its maximum near the wall, attaining there a value 
on the order of 55% and decreasing very rapidly with increasing 
distance since the f^ frequency was dominant. 
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Figure 10 

Ordinate Dependency of 
Energy of Instabilities at 
Frequencies f]^ (a) and 
f 2 (b) and of the Har- 
monic Rate (c) 

U„=16 m/sec, 2a=-l mm, 
x*=50 mm, 63 = 2 . 27 


We were especially interested 
in the maximum energy level for fre- 
quency f 2 and in the maximum value of 
the energy ratio between the two fre- 
quencies . 

6.2. Variation as a Function of the 
Abscissa 


The set of curves in figure 11 
represents the dependence of max- 
imum energy, i.e. the quantity 

(/u ' 2 /u ) ^ , on X* for various 

e t2max 

deformation amplitudes. Notice that 
the greater the amplitude of the 
deformation, the more visible the 
phenomenon was. Note also that the /435 

curves have a maximum at an abscissa 
X* of between 50 and 70 mm, depending 
on the deformation amplitude, and 
therefore downstream from the middle 
of the cavity. 



:'2 vs. Abscissa for Various Wall {Uoo = 16 m/sec) 

Deformation Amplitudes (Uoo=16m/sec) 
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On the contrary, curves representing the variation of the 
maximum ratio between the two energies as a function of x* (fig- 
ure 12) attain a maximum for abscissas of between 40 and 60 mm, 
and therefore upstream from the middle of the cavity. The curves 
then rapidly tend toward horizontal asymptotes. 



'■Max. ■ 

V. 

30 


20 


to 


Figure 13 

Dependence of Energy of 
Instabilities at f 2 (a) 
and of Harmonic Rate (b) 
on Wall Deformation 
Amplitude (2a) at x*=45 mm 
Uoo = 16 m/sec 


In figure 13, we have 
plotted the curves showing the 
dependence on amplitude 2a of max- 
imum energy at frequency f 2 (curve 
a) and of the maximum energy 
ratio between the two frequencies 
(curve b) . The latter curve exhibits 
a change in curvature between 
-0.75 and -1 mm, which does not exist 
on curve a. This clearly shows the 
difference in behavior between 
frequency f 2 and the ^^"^2 
ratio . 


6.3. Discussion of the Harmonic's Origin 


On the basis of these experimental results, various proposals 
could be advanced to explain the harmonic's origin. 

We have already eliminated the hypothesis of a distortion 
arising in the processing of the hot wire signal, but it remains 
possible that the physical phenomenon itself is nonlinear. This 
would result in the propogation of harmonics. 

However, this proposition does not seem plausible when the 
harmonic ratio attains values as high as 50%. Furthermore, let us go 
back to figure 12 and consider for example the variation in the 
harmonic ratio between x* = 50 mm and 60 mm for a -1 mm wall defor- 
mation. We find the following values: 
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X* == 50 mm 


X* = 60 mm 


(Vt/VU,)/. 4.96,10-5 22.4.10-5 

2.6.10-5 3.6.10-5 

(%//7'VU,)V('//rVU,)/. 52% 16,2% 

These values show that while the fundamental frequency is 
amplified by a factor of about 4, the harmonic ratio diminishes 
by 2/3. This is not compatible with the hypothesis linking the 
presence of the instabilities' first harmonic to the nonlinear 
character of their development. If a nonlinear phenomenon were 
in fact involved, the harmonic ratio would increase in a con- 
tinuous manner when the fundamental frequency was amplified. 

If in contrast we look at the geographic location of the region 
where the harmonic existed, we observe that the zone where the har- 
monic ratio was greatest corresponds to the zone where the exterior 
velocity gradient was the most negative. 

For a zero gradient, we have seen (figure 6) that the f^^ 
frequency corresponded to the frequency having m.aximum amplifica- 
tion according to stability theory. Now, a negative gradient 
provoked a widening of the range of unstable frequencies and if 
the gradient was sufficiently negative, the £2 frequency fell in 
the unstable area. 

To give an example, for the abscissa x* = 50 mm, corresponding 

in a flat plate (A = 0) to ^ = 850, the frequency with maximum 

amplification had a reduced pulsatance w* of 1.1 x 10 In 

contrast, in the conditions of flow (A = -12) the pulsatance 

-4 

corresponding to the harmonic, oi* = 2 x 10 , had maximum amplification 

and the fundamental was located in the unstable area. The ratio 
of about 2 between the two pulsatances led us to make the hypothesis 
that the harmonic encountered in the experiments was due to a 
laminar microflow along the deformation. 

We then used the split film sensor to measure the velocity 
vector's angles of incidence in the boundary layer. Figure 14 
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depicts the velocity vector along 
the deformation (magnitude and 
direction) . There did not appear 
to be any return flow in the 
neighborhood of the wall. There 
were, of course, strongly per- 
turbed zones, but the velocity vec- 
tor's angle of incidence always 
remained relatively small. 

Under these conditions, we 
rejected the hypothesis of a permanent flow as the origin of the £2 
frequency. However, observation of the instantaneous signal showed 
us the intermittent character of ^ 2 '^ appearance in the unstable 
bursts. This led us to readjust our hypothesis by proposing that 
the f 2 frequency could be due to an intermittent flow [10] . The 
absence of a threshold in the variation of the harmonic and the 
harmonic ratio with deformation amplitude (figure 13) could then 
be explained by a gradual increase in the intermittence of the 
microflow. 

7. Conclusion 


V 

(mm) 

1 . 

.8 

.6 

.4 

.2 

0 


Figure 14 

Representation of Velocity 
Vector in the Wall Deforma- 
tion Zone, 2a = -1 mm 
(U<„ = 16 m/sec) 



The natural instabilities which arise in a flat plate's 
laminar boundary layer are composed of intermittent wave trains. 
Spectral analysis of these fluctuations localized their frequencies 
and allowed us to isolate their amplitude modulation spectrum. 

Although there is a fundamental difference between natural insta- 
bilities and those under consideration in Tollmien's and Schlichting ' s 
theoretical calculations , the frequencies of the instabilities and 
their amplification coefficients were identical to those predicted 
by stability theory. A wall deformation does not change the 
structure of the phenomenon. It only induces an exterior velocity 
gradient which, according to its sign, either amplifies or attenuates 
the instabilities to a degree depending on their energy level. Obser- 
vation of the instantaneous velocity fluctuation signal as well as 
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spectral analysis showed that in a zone of negative gradient the 
first harmonic of the instabilities appears in the form of a perturba- 
tion grafted onto the unstable bursts of frequency f^. 

Given the significance of the harmonic under certain experimental 
conditions, its existence could not be satisfactorily explained by 
a nonlinear characterization of the phenomenon. 

We also eliminated the hypothesis of a permanent microflow being 
at the origin of the perturbation at frequency £ 2 - In contrast, 
the intermittent character of this perturbation is an indication 
that it might be due to a microflow that is itself intermittent. 

Such a phenomenon could not be detected by the split film anemom- 
eter sensor, which only provided average values for the velocity's 
angle of incidence. 

Apart from the instantaneous signals and their spectra, all 
the measurements we have presented correspond to temporal averages. 
This gives the study a generalized character. The problematic 
intermittence of the instabilities makes a statistical analysis of 
the phenomenon necessary, in fact. We will thus attempt to 
restructure the boundary layer by imposing a definite, known 
artificial intermittence on it. To do this, we are investigating 
the influence of sinusoidal wall vibrations on the boundary layer. 

(Manuscript received October 30, 1975, revised February 22, 1977.) 
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